Experimental evidence of accelerated energy transfer in turbulence 
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We investigate the vorticity dynamics in a turbulent vortex using scattering of acoustic waves. 
Two ultrasonic beams are adjusted to probe simultaneously two spatial scales in a given volume of 
the flow, thus allowing a dual channel recording of the dynamics of coherent vorticity structures. 
Our results show that this allows to measure the average energy transfer time between different 
spatial length scales, and that such transfer goes faster at smaller scales. 

PACS numbers: 47.32.Cc 43.58. +z 47.27.Jv 



Much of the investigation in turbulence has been de- 
voted to the understanding of the mechanisms underly- 
ing the energy transfer through the turbulent cascade, a 
concept introduced by Richardson in 1922 [l|. A good 
deal of theoretical and numerical work has been dedi- 
cated to study shell models, which include the essential 
features expected in a turbulent cascade, but without 
some of the complications inherent to the Navier-Stokes 
equations. Though these models retain some of the dy- 
namics of the motion equations, they handle the velocity 
field as a scalar. Thus the price payed is a complete loss 
of flow geometry Q. Another approach is the study of 
random multiplicative cascade models, which adequately 
mimic the statistics of the local energy dissipation rate 
e at different flow scales, and the intermittency of small 
scales observed in a large number of experimental results. 
Although much progress was done, many issues related 
to the statistics of energy transfer, energy dissipation and 
small scale intermittency in turbulent flows still remain 
without satisfactory answers. In particular, the dynamics 
of the energy transfer through the turbulent cascade has 
been only recently addressed in some theoretical works 
[H 0); but this aspect of turbulent flows is still lacking 
experimental studies specifically addressing it. 

In this note we report an experimental result, obtained 
for the case of a vortex embedded in a turbulent flow pro- 
duced in air by two coaxial centrifugal fans, each facing 
the other and rotating at fixed angular velocities. This 
so-called von Karman flow has been used in a number of 
experiments in turbulence [5(. When the fans rotate in 
the same direction with identical speeds, a strong vortex 
is produced between them Q . The advantage of this con- 
figuration is that in a small volume we obtain an intense 
and well-defined coherent vorticity structure, surrounded 
by a fully turbulent background. Details of this setup are 
given in a previous work [7| . The parameters used in this 
case are: diameter of fans: D = 30 cm, height of vanes: 
h = 2.2 cm, distance between disks: H = 30 cm, rota- 
tion speed: / = 30 Hz, and diameter of central holes: 
d = 2.5 cm. The fans were driven by two DC motors, 
powered by independent constant voltage sources allow- 
ing to keep the disk rotation within 1% of the desired 



speed. Local measurements of the airflow speed were 
performed with a Dantec Streamline 90N10 frame hous- 
ing a Constant Temperature Anemometer which drives a 
5 /im diameter and 2 mm length platinum wire probe. A 
Streamline 90H02 Flow Unit was used to calibrate it. To 
produce and detect the sound waves, four circular Sell 
type electrostatic transducers [|[ were used, each having 
an active surface diameter of 14 cm. Two B&K model 
2713 power amplifiers drove the emitters. The receivers 
used custom made charge amplifiers, whose output were 
connected to two Stanford Research model SR830 lock-in 
amplifiers (LIA) through 2 nd order high-pass RC filters 
having cut-off frequencies of 7.5 kHz. The sine-wave out- 
put of each LIA internal generator was used as signal 
source for the corresponding emitter. The LIA output 
signals, comprising in-phase and quadrature components, 
were low- pass filtered by an IOTech Filter-488 8 th order 
elliptic filter bank, using a cutoff frequency of 4 kHz. The 
resulting signal was digitized at 12.5 kS/s with a 16-bit 
resolution National Instruments AT-MIO-16X multifunc- 
tion board, installed into a personal computer. A total 
of 2 20 points per signal channel was collected in each of 
77 data sets. 

Fig. 1 (top) shows the device that produced the flow, 
along with the four circular ultrasonic transducers, and 
the hot-wire probe. The transducers were arranged 
to perform measurements using an interferometric con- 
figuration, described in detail in a previous work [8[. 
The scattering angles were Q\ = 40° and 62 = 60°. 
The magnitude of the scattering wave-vectors, q? catt = 
47w , 4 sin(6>f at 72)/c = \k' - k\,i = 1,2, were adjusted 
independently for each channel through their correspond- 
ing frequencies vq^. We acquired three sets of data, set- 
ting 1/0,1 = 20 kHz, 40 kHz, and 60 kHz; while i/o,2 took 
some 25 values in the interval [16, 90] kHz in each case. 
Notice that each probed scattering wave-vector, related 
to the z -component of the vorticity distribution, can be 
set arbitrarily close (even identical) to the other. Fig. 1 
(bottom) displays a schematic upper view of the lower 
disk, hot-wire probe and ultrasonic transducers, along 
with distances and angles. This arrangement is symmet- 
ric with respect to the vertical plane A A. Also displayed 




FIG. 1: Top: Partial view of the experimental setup. Bottom: 
a schematic drawing of the experimental setup (see text). 



is a schematic diagram of the measuring chain. Fig. 2 
(top) is a sketch of the volume probed by the ultrasonic 
beams, showing the relationships between scattering vec- 
tors and wave-vectors in the case of exact inter-channel 
tuning (same scattering wave- vector) . As mentioned in 
previous works @, 0], the vortex performs a slow pre- 
cession motion around the axis of rotation of the fans. 
Even so, it remains inside the measurement volume, as 
shown in fig. 2. A horizontal cut of the measurement 
volume is displayed as a shaded rhombus, the heavier 
shaded circle indicating the position of the vortex core. 
Air speed measurements with the probe wire aligned ver- 
tically, and located at about 2 cm of the disks's rotation 
axis and some 10 cm above the lower disk, give the sam- 
ple record displayed in fig. 2 (bottom). The nearly peri- 
odic precession of the vortex is visible as a pattern with 
sharp dips — sometimes enlarged or doubled by fluctu- 
ations due to the turbulent background, each time the 
center of the vortex is close or coincides with the hot 
wire location. These events are signaled by arrows in 
fig. 2, the time interval between single events being about 
1 s. Accordingly, the spectrum displayed in fig. 3 (top) 
has a rather wide peak located at /o ~ 1 Hz. The Tay- 
lor hypothesis is not met by this flow, but this spec- 
trum gives an idea of the energy contents at different 
scales. At higher frequencies (smaller scales), it meets 
only approximately the fc~ 5 / 3 K41 law (represented by 
the straight line) in the inertial range, due to the pres- 
ence of the vortex in the bulk of the flow. In fact, we 
observe three regions, A, B, and C, where the spectrum 
looks self-similar, but with different slopes. As we will 
see later, this could be related to the influence of the 
flow anisotropy on the dynamics of the energy transfer 




FIG. 2: Top: Diagram showing the relationship between 
wave- vectors, the volume probed by the ultrasonic beams 
(shaded region) and the size of the disks. The dark circle 
represents a cut of the vortex core. Bottom: A sample of 
the anemometer signal. Arrows signal events of minimal air 
speed, the vortex core being close to the hot wire probe there. 
The inset displays the velocity profile of a non turbulent vor- 
tex, which helps in recognizing the vortex core in the turbu- 
lent signal. 



through the turbulent cascade. The acoustic signals are 
delivered by the LIAs as low frequency complex voltages, 
which are images of the spatial Fourier modes U, z (qj , t) 
of the flow vorticity, probed at well-defined wave vectors 
qj = 47ri/ j sin(#j/2)/c, j — 1, 2. Typical spectra of these 
signals are displayed in fig. 3 (bottom). The wide central 
peak, due to the intense coherent vorticity of the vortex, 
lies between two side bands produced by the background 
flow. There is a close correspondence between the low fre- 
quency region of the velocity spectrum (up to 50 Hz) and 
the central region of the vorticity spectrum. The same 
is true for the high frequency region (> 200 Hz) of the 
spectra, both exhibiting a slowly decaying roll-off, with a 
slope close to —5/3 related to the background turbulence 
of the flow. 

Now we address the separation of the coherent 
structure and the background turbulent flow. Fol- 
lowing ref. Q, we compute the coherence between 
two spatial Fourier modes of the vorticity, Q, z (qi , i) 
and £l z (q2,t). The coherence function, C(q±, q2, v) = 
(tt z (qi,v)n* z (q 2 , z/)| 2 ><|fi*(g 2 , z/)| 2 >, where 

v is a frequency shift associated to the Fourier transform 
in time of tt z (qi , i) , is a statistical estimator that allows 
to evidence phase stationarity (coherence) between 
distinct Fourier modes. Fig. 4 (top) is a plot of the co- 
herence between two scattering channels tuned to probe 
the same scale of the vorticity distribution (q\ = q^). A 
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FIG. 3: Top: Power spectral density of the signal displayed 
in fig. 2. Three almost linear regions are labeled A, B, and 
C (see text), with slopes q.a,b,c- The straight line represents 
the Kolmogorov spectrum, with slope olk = —1.66. Bottom: 
Individual spectra of the signals delivered by the lock-in am- 
plifiers. The spectrum corresponding to the scattering angle 
of 60 degrees decreases slightly faster with the frequency shift. 
Scattered waves having the same central frequency of the in- 
cident beam have no frequency shift. 



central peak with a high level of coherence (0.9, being 1 
the maximum expected value; the value of 0.9 should be 
ascribed to incoherent noise) reveals that both channels 
are detecting the same coherent structure. This was 
verified "turning off" the vortex by rotating the disks 
in opposite directions at the same angular speed. In 
this case only a turbulent shear layer remains, and the 
plot of the coherence function displays only the side 
bands. By stacking plots of coherence for a sequence of 
values of q\, it is possible to build surfaces of coherence 
relating events at different flow scales, as functions of 
Aq = <72 — qi and the frequency shift v, as shown in fig. 4. 
We can interpret them in as eddies at given scales draw- 
ing its energy from decaying structures at larger scales, 
in the frame of the Richardson's turbulent cascade. As 
these surfaces represent a time average of the whole 
story, we perform a time cross-correlation analysis to 
confirm this interpretation. This is done by computing 
the cross-correlation x(<?i, Aq, At) = (\il z (qi,t)\\il z (qi + 

A q ,t + At)\}/^(\n z ( qi ,t)\2)(\n z ( qi + Aq,t + At)P>, 

were \ fl z (qi, t)\ is a measure of the strength of a structure 
with vorticity £l z (qi,t). The function Aq, Ai) will 
have a peak at At = r (r = time lag), indicating the 
presence of a vorticity pattern at the scale l\ = 2n/qi 
in time t that appears at the scale Z 2 = lnjqi a time 
t + t later. The peak amplitude is a measure of how 
well the pattern is preserved from one scale to the other, 
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FIG. 4: Top: coherence plot of signals detected simultane- 
ously by both receivers. The wide peak centered at zero fre- 
quency shift is due to the vortex. Bottom: Three coherence 
surfaces are shown, obtained by detuning one of the emit- 
ters by a frequency corresponding to a shift Aq in the wave 
number. 



a value of one meaning a perfect preservation. Fig. 5 
shows three straight lines resulting from fits to the time 
lags calculated from such cross-correlation analysis, each 
corresponding to one of the coherence surfaces shown 
in fig. 4. The time lags found here are consistent with 
those given in reference Q. We have r ~ 1 ms for 
Aq = 8 cm -1 , and the transfer rates of fig. 7 (c) in 
reference Q are similar [loj . 

Our main result is the experimental observation of ac- 
celeration in the energy transfer in turbulence, evidenced 
by the systematic change in the trend of the data in fig. 
5 when the magnitude of the reference wave-vector q\ 
is changed. The slopes of the fitted straight lines de- 
crease with increasing q\. The scattering of the points 
is an indication of the intrinsic variability of the energy 
transfer time through the turbulent cascade (R 2 > 0.6 in 
all cases). The following argument helps to understand 
the trend change: the time lag measured between the 
wave numbers q\ and (72 = 9i + Aq can be seen as the 
average transfer time of the energy contained in a scale 
corresponding to q\ to a smaller scale corresponding to 
q-2- Thus, the slope Sr/Sq (we drop A's from now on) 
of the straight lines can be related to the inverse of the 
average transfer rate of energy in the k -space, in a neigh- 
borhood of the reference scale corresponding to the wave 
number q\ [11]. This decrease in slope with increasing 
qi reveals that the average transfer time of the energy 
decreases when the involved scales are smaller. In other 
words, taking 5q/Sr as the average transit speed (in the 
k -space) of structures starting at the wave number qi 
and ending at the wave number q2, we conclude that 
such speed is larger at larger q's (or smaller scales). This 
is consistent with a remark made in a previous numer- 
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FIG. 5: Time lags corresponding to the coherence sur- 
faces of fig. 4. a), b), and c) correspond to qi = 
3.79, 7.58, and 11.4 cm -1 , respectively. The straight lines 
are linear fits to the measured points, with correlation coef- 
ficient R 2 = 0.60 (0.75,0.66). Note that the slopes decrease 
with increasing q\. 

ical work [3(. A rather crude explanation of this result 
follows: if E{q) is the spectral density of energy per unit 
mass — averaged in time, and we think of Sq/Sr as be- 
ing an ordinary derivative, then E(q)5q/Sr is the average 
energy flux density in k -space. But this is the average 
dissipation rate per unit mass e. Thus, E(q)5q/ 'Sr = e, 
and Sr/Sq — E(q)/e. As E(q) decreases with increasing 
q, so does Sr/Sq. Notice that in this picture, the energy 
attached to a given eddy is transported towards larger 
wave numbers merely by reduction of its spatial scale (in 
directions perpendicular to its vorticity) due to stretch- 
ing by the velocity field of larger eddies. 

We can illustrate our view in the frame of K41 the- 
ory [l|: we have E(q) = Cke 2 / 3 ^ 5 / 3 , where Ck is 
the Kolmogorov constant. Thus, for homogeneous and 
isotropic turbulence, Sr/Sq — Cke -1 / 3 ^ 5 / 3 , which is of 
course a decreasing quantity. By rearranging terms, we 
have q~ 5 ^ 3 Sq = (e 1 / 3 /Ck)6% which by integration yields 

q(t) = (q 2/3 - fe^t/Ckr 3 / 2 . This is the itinerary of 
the average magnitude q in a homogeneous and isotropic 
turbulent flow in which the mean dissipation rate per unit 



mass is e. The scale Iq = 2ir/qo, corresponding to t = 0, 
can take any value within the inertial range. Notice that 

q{t) diverges in a finite time t K — ^Ck e^^q^ 2 ^ , corre- 
sponding to the average time taken by the energy at the 
scale Zo to cascade down to I = + (q — > oo), in a fluid 
with vanishing viscosity — in the frame of K41 theory. 
Putting Ck ~ 0.5, e ~ 10 2 W/kg, estimated from the 
power injected to the flow (« 70 W) and the volume of 
moving air (~ 0.5 m 3 ) and, for instance, qo ~ 3 x 10 2 
m _1 , corresponding to the largest scale measured with 
our acoustic method, we get t K ~ 4 ms. This is on 
the order of characteristic times in fig. 5, showing that 
transfer times in homogeneous turbulence do not differ 
too much from those in our flow. Thus, our reasoning 
shows that the trend observed in fig. 5 can be understood 
in terms of "classical" homogeneous and isotropic tur- 
bulence. In non-homogenous, non-isotropic turbulence, 
as in the present experiment, the behavior of Sr/Sq can 
be estimated from the experimental spectrum shown in 
fig. 3, where the Kolmogorov spectrum is plotted as a 
straight line, for comparison. In the large scales region 
(A), Sr/Sq decreases faster than in K41 turbulence, the 
cascade acceleration being larger. In the small scales re- 
gion (C), Sr/Sq tends to reach the Kolmogorov slope, 
but the acceleration is slightly smaller than in K41. Re- 
gion (B), in which the cascade acceleration reaches its 
minimum, is a transition zone from a non-isotropic flow 
at large scales, where vorticity is mainly parallel to the 
z-axis, to a more isotropic small-scale flow. In all three 
regions the spectrum seems to be self-similar, with dif- 
ferent scaling exponents. Of course, being this one a 
first experimental result revealing the accelerated nature 
of energy transfer through the turbulent cascade, more 
work will be necessary to state our conclusions on more 
solid ground. We recall that we are measuring only the 
z -component of the vorticity and the modulus of the ve- 
locity component normal to the z-axis. It should be 
necessary to simultaneously measure at least one more 
component of the vorticity. Additionally, the conditions 
of validity for the Taylor hypothesis are not meet in these 
hot-wire measurements — a flying probe would be much 
better to study this flow. At any rate, in the large-scale 
range the vorticity in our flow is mainly parallel to the z - 
axis. Thus, eddies in the large scale region have mostly 
parallel vorticities, and the vortex stretching is largely 
inhibited there, slowing down the energy transfer. At 
smaller scales, where a more isotropic flow exists, vortex 
stretching is reestablished, giving faster energy flux rates. 
Thus, in-between the acceleration should be greater than 
in the isotropic case. This could explain in part our re- 
sults. To conclude, we remark that our reasoning suggest 
that cascade acceleration may be a necessary but not suf- 
ficient condition for small scale intermittency. 
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